Importance of local correlations for the order parameter of high-T c superconductors 
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It is shown that local temporal correlations in addition to non-local spatial correlations are im- 
portant to understand the size and the doping dependence of the d-wave superconducting order 
parameter of high-temperature superconductors. To this end, the hole- and electron-doped two- 
dimensional Hubbard model at zero temperature is considered and treated by an extension of the 
variational cluster approximation. Within this approach, the effects of temporal correlations can be 
studied systematically and in a thermodynamically consistent way by comparing results obtained 
from different reference clusters. Contact can be made with previous cellular (plaquette) dynami- 
cal mean-field calculations. This shows that temporal correlations considerably decrease the order 
parameter and provide a substantial gain of binding energy. Besides, a few methodical insights 
regarding real-space quantum-cluster approaches are obtained in addition. 

PACS numbers: 71.10.-w, 74.20.-z 



I. INTRODUCTION 

There are strong efforts devoted to a convincing the- 
ory of high-temperature superconductivity with non- 
local d-wave order parameter. Some agreement has been 
achieved that essential physical properties can be stud- 
ied within the prototypical two-dimensional single-band 
Hubbard model. 1 Using standard notations the (grand- 
canonical) Hamiltonian reads: 



H = 



(Uj -M%)c| ff c iCT + U 



(1) 



with c\ a (da) being the creation (annihilation) opera- 
tor of an electron at lattice site with spin projection 
a =t)4- n ia = c\ a Ci a denotes the corresponding occupa- 
tion number operator. The energy scale is set by the 
nearest-neighbor hopping t nn = — 1. To describe the 
universal low-energy physics of cuprates, an additional 
second-nearest-neighbor hopping term with t nnn = 0.3 2 
as well as U — 8 are taken as standard parameters which 
will be kept fixed throughout the paper. 

The physics of the intermediate- to strong-coupling 
regime of the model is notoriously complicated since 
there are several phases with different long- or short- 
range correlations competing on a low-energy scale. Due 
to the non-locality of the (d-wave) order parameter in 
the superconducting phase, a pure (single-site) mean- 
field approach cannot capture the essential features of 
the model. On the other hand, a direct numerical 
solution of finite but large two-dimensional Hubbard 
lattices is in principle possible via a quantum Monte- 
Carlo approach. 3,4 However, for the interesting parame- 
ter regime, i.e. doped systems at low temperatures, the 
so-called sign problem makes simulations ineffective. 

In this situation, quantum-cluster approaches^ ap- 
pear promising, i.e. cluster extensions^— of the dy- 
namical mean- field theory (DMFT),— ^ or the varia- 
tional cluster approximation (VCA) i 12 ' 13 Common to 
all quantum-cluster approaches is that, for the calcu- 



lation of the electronic self-energy, a small cluster is 
self-consistently or variationally embedded in a non- 
interacting bath that approximately accounts for the ef- 
fects of the cluster environment. 

Two types of correlations must be distinguished here: 
(i) Spatial correlations are neglected completely in single- 
site DMFT but are included within a cluster approach, 
up to the size of the cluster. Here the number L c of cor- 
related sites with finite Hubbard- U in the cluster is es- 
sential. For the doped two-dimensional Hubbard model, 
the important feedback of non-local magnetic correla- 
tions on the single-particle spectrum, for example, can 
only be captured by a cluster approach. As spatial cor- 
relations are neglected beyond the size of the reference 
cluster, a quantum-cluster approach can also be seen as a 
cluster mean-field theory, (ii) Temporal correlations are 
fully accounted for already in single-site DMFT and give 
rise to a highly non-trivial frequency dependence of the 
self-energy which, for example, is vital to understand the 
Mott transition. Here the number of uncorrelated sites 
with U — representing the bath is essential. As in the 
DMFT a continuum of bath degrees of freedom is used, it 
can be considered as an optimal mean-field theoryfiSr— 
If a continuum of bath sites is used, as it is intended 
usually, 7-9 an optimal quantum-cluster approach is gen- 
erated for a given L c . 

Besides quantum Monte-Carlo techniques, the Lanc- 
zos method^ is frequently used as a cluster solver since 
this allows to study the ground-state phase diagram zero 
at temperature T = 0. As the Lanczos method is lim- 
ited by the total number of sites, the variational clus- 
ter approximation (VCA), if combined with Lanczos as 
a solver, usually takes into account as many correlated 
sites as possible but completely disregards bath degrees 
of freedom i 12 ' 13 The idea is that for a large cluster, tem- 
poral correlations are sufficiently accounted for since they 
are restored anyway in the infinite-cluster limit L c — > oo 
where the VCA (as any quantum-cluster approach) be- 
comes exact. 

Using the VCA (no bath sites included), the ground- 
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state phase diagram of the two-dimensional Hubbard 
model has been explored using clusters with up to L c = 
12 sites and different cluster geometries.— ~— Close to 
half-filling, for hole as well as for electron doping, these 
calculations reveal a phase where antiferromagnetic order 
and d-wave superconductivity are microscopically coex- 
isting (AF+SC). With increasing doping a pure super- 
conducting phase (SC) persists in both cases. In agree- 
ment with experimental data, the phase diagram is asym- 
metric: Antiferromagnetic order, for example, extends to 
higher doping in the electron-doped case. For the hole- 
doped system, an extended macroscopic phase separation 
of the mixed AF+SC phase with a purely superconduting 
phase SC can be found.— This has been interpreted as a 
tendency towards the formation of microscopically inho- 
mogeneous (e.g. stripe) phases. On the other hand, there 
is hardly a phase-separated state for the case of electron 
doping which may be seen to be consistent with the ab- 
sence of rigorous signs for stripe structures in electron- 
doped materialsi 16 i 18 

Using the cellular DMFT (including bath sites) for the 
set of hopping parameters considered here^S but also 
for the particle-hole symmetric case (t nnn = 0)j2! one 
can nicely reproduce the dome-like shape of the super- 
conducting order parameter found in experiments. For 
tnnn = 0.3 a homogeneous coexistence of antiferromag- 
netism and superconductivity (AF+SC) is found for fill- 
ings close to half-filling, while the pure superconducting 
phase (SC) persists for larger hole and electron concen- 
trations. For the hole-doped system the pure supercon- 
ducting phase extends over a broader range of dopings 
than on the electron-doped side. 

We conclude that VCA (no bath sites) and C-DMFT 
(including bath sites) yield very similar ground-state 
phase diagrams. However, there is an obvious discrep- 
ancy with respect to the size of the superconducting order 
parameter. VCA and C-DMFT results can differ by more 
than a factor 2 (see also results below). In addition, the 
optimal (hole and electron) doping is larger in the VCA 
as compared to the C-DMFT calculations. A major pur- 
pose of the present paper is to point out that this is due to 
an underestimation of local temporal correlations within 
conventional VCA. Physically, this means that temporal 
correlations, giving rise e.g. to Kondo screening of mag- 
netic moments, are important to understand the order 
parameter, the phase diagram and eventually the pairing 
mechanism. 

For practical calculations using the Lanczos method it 
is quite tempting to consider a plaquette of four corre- 
lated sites since this allows for d-wave order as well as 
non-local singlet formation with a minimum computa- 
tional effort. Employing a plaquette represents the first 
important step beyond a single-site (dynamical) mean- 
field approach. Here we present calculations obtained by 
an extension of the conventional VCA, i.e. we employ a 
reference cluster with one additional bath site attached to 
each correlated site. This can be done within the frame- 
work of the SFT by comparing different but thermody- 



namically consistent approximations. 

II. SELF-ENERGY-FUNCTIONAL THEORY 

Both, the VCA and the C-DMFT, can be seen as 
approximations originating from a certain cluster refer- 
ence system within the context set by the self-energy- 
functional theory (SFT) ^^2, I n the (conventional) VCA 
we consider a finite cluster without bath sites embed- 
ded in the original lattice while for C-DMFT a continu- 
ous bath is attached to each correlated site of a cluster 
(see Fig. [I}. Both include non-local short-range spatial 
correlations on a scale up to the linear size of the clus- 
ter. Beyond that both approximations are mean-field- 
like. The treatment of local (temporal) correlations, how- 
ever, is different: The continuous bath considered in the 
C-DMFT ensures that these are taken into account ex- 
actly for arbitrary cluster size and even in the case of 
single-site DMFT. On the other hand, within the VCA 
the description of local correlations becomes exact in the 
infinite cluster limit only. 

Self-energy-functional theory^ starts from the grand 
potential of a system of interacting electrons expressed 
as a functional of the self-energy, 

i2[S] = Trln(G c 7 1 -S) _1 + F[S] , (2) 

with the free Green's function Go and F[E] being the 
Legendre transform of the Luttinger-Ward functional 
<f>[G].— This functional can be shown to be stationary 
at the exact (physical) self-energy. Hence, we have the 
following dynamical variational principle: 

<Jfi[£] = 0. (3) 

This, however, cannot be evaluated in practice since the 
functional form of F[S] (and of $[G]) is actually un- 
known. The main idea of the SFT is to restrict the vari- 
ation of the self-energy in the variational principle Eq. 
@ to a certain subspace of trial self-energies which is 
spanned by the self-energies of an exactly solvable ref- 
erence system (i.e. a small cluster). This means to pa- 
rameterize the trial self-energy Sf by the one-particle 
parameters of the cluster t' and to treat t' as variational 
parameters: 

I^Pt'] - 0. (4) 

For a small cluster, the value of the grand potential at 
the trial self-energy, f2 [£$/], and thus the condition Eq. 
(|4]) can be evaluated numerically exact, see Refs. [l2lfl8l 
for examples. 

For any practical calculation, one usually considers a 
few physically relevant parameters only. Fig. Q] illustrates 
the different plaquette reference systems considered in 
this study and, for the VCA, the one-particle parame- 
ters to be optimized. The conventional VCA refers to a 
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VCA C-DMFT 

□:• ZZ 

continuous bath 




VCA + bath sites 

FIG. 1: Reference systems and variational parameters as used 
for our calculations. Additionally a corresponding reference 
system generating the C-DMFT is illustrated. Here, a contin- 
uous bath is optimized while the parameters of the correlated 
sites remain at their physical values. 



reference system with the following Hamiltonian: 

H> = 53 t'iA* ^ + ( £ ~ A*) 53 + 17 53 

ija ia i 

+ h sc 53 A»j (cnc jt + H.c.) (5) 

ij 

including a global shift of the on-site energies e and a fic- 
ticious symmetry-breaking Weiss field of strength hsc as 
variational parameters. Optimization of e ensures ther- 
modynamic consistency with respect to the total particle 
number 17 while the Weiss field allows for a possible super- 
conducting phase. Note that the interaction term is the 
same as for the original model Eq. (JTJ). Contrary, only 
intra-cluster hopping parameters iy are retained while 
the inter-cluster hopping is switched off. For d x i_ y i- 
pairing, 

A = / 1 for n~r j = ±e x 
13 y -1 for r t - Tj = ±e y v ' 

with cluster sites r^, Vj and e Xl e y being the unity vectors 
in the x and y direction. 

Attaching one additional bath site to each correlated 
site results in the reference system Fig. [TJ bottom. The 



Hamiltonian reads: 

H ' = XX^ia^ + ( £ - M)X] ni °- +t7 X] rii t ri ^ 

ija ia i 

+ V 53 ( a l<r c ^ + H.c.) + (e b - m) 53 nh ^ 

ia ia 

+ h h ,SC 53 ( a il a Jt + H -°-) • ( 7 ) 

ij 

Here, the hybridization strength V with the bath sites 
and their on-site energy Eb are treated as additional vari- 
ational parameters. This leads to an improved descrip- 
tion of local correlations as compared to the reference 
system Eq. J5]). In Eq. a\ a (a ia ) creates (annihilates) 
an electron with spin projection a =f , | on a bath site 
coupled via V to the cluster site r,, and nh,ia = a itJ ai a . 

Within the SFT and using the reference system in Fig. 
UJ (right), it turns out that the Euler equation of the vari- 
ational principle Eq. (QJ is equivalent with the C-DMFT 
self-consistency equation for the parameters of the con- 
tinuous bath. The optimal values of the one-particle pa- 
rameters associated with correlated sites can be shown to 
be given by their physical values, i.e. those of the original 
systemi^ 7 - For a reference system with a finite number of 
bath sites, on the other hand, this is no longer necessarily 
true. Namely, the truncation of the inter-cluster hopping 
is partially compensated for by optimal parameters as- 
sociated with the correlated sites that differ from their 
physical values. In the conventional VCA (i.e. no bath 
sites) this is obvious. 

For a refererence system including bath sites, Fig. [JJ 
(bottom) , we therefore expect that applying the ficticious 
(i.e. unphysical) Weiss field to the bath sites is much more 
efficient as compared to the correlated sites. We have nu- 
merically checked this for the case of antiferromagnetic 
order where both, a staggered magnetic Weiss field on the 
correlated and on the bath sites have been varied simul- 
taneously: The optimal staggered field on the correlated 
sites has turned out to be typically more than an order 
of magnitude smaller than the one on the bath sites and 
to be negligible for the calculation of observables. 

Consequently, in Eq. we attach the symmetry 
breaking Weiss field /ib,sc to the bath sites. Note, how- 
ever, that we still do optimize the parameter e (in ad- 
dition to Eb) to ensure the above-mentioned consistency 
with respect to the particle number. 

III. SPIN-DEPENDENT PARTICLE-HOLE 
TRANSFORMATION 

As has been demonstrated by Senechal et alJ^ a 
possible way to treat superconductivity within VCA is 
to employ the Nambu formalism where normal as well 
as anomalous Green's functions have to be computed 
to evaluate the self-energy functional. An alternative 
is given by a spin-dependent particle-hole transforma- 
tion of the original and the reference system. This re- 



4 



stores particle-number conservation and avoids anoma- 
lous Green's functions. The transformation is given by 



J 



4r 



n 



(8) 



Applying the particle- hole transformation to H' , Eq. (J7J, 

we get 



H' = 



+ (e - M) X! ("*t - "ij.) + 17 X! 

i 

+ ^ (e + e b - 2fi) , 



(9) 



where ru 



did,., 



have been used. 



6j b ia and ^, 



'" ~~ ^ aild 

Ay = Aji 

Eq. ([9]) is interpreted as a reference system correspond- 
ing to the particle-hole transformed original system Eq. 
((T|). For both cases, the transformation yields an attrac- 
tive Hubbard interaction (U — > —U). Hopping terms be- 
come spin-dependent, and chcmical-potcntial-like terms 
are mapped onto ferromagnetic fields and vice versa. The 
nonlocal pairing field maps onto a nonlocal spin-flip term. 
The transformed Hamiltonian respects particle-number 
conservation (unless there were spin-flip terms in the 
original one which would transform to pairing fields), 
the z-component of the total spin, however, is no longer 
conserved. For the transformed reference system, the 
anomalous Green's function vanishes; the normal Green's 
function is no longer diagonal in the spin index a. 

In practice, we first evaluate the self-energy functional 
for the transformed original system H by solving the 
problem posed by the transformed reference system H', 
i.e. we calculate the corresponding Green's function and 
self-energy, evaluate therewith the self-energy functional 
and optimize the variational parameters. Note that the 
variational parameters are chosen accordingly (e.g. the 
on-site energy is represented as the strength of a ferro- 
magnetic field in the transformed system, see Eq. ^j, 
and this field strength is varied and optimized). We have 
crosschecked with the results obtained without transfor- 
mation where possible. Observables, i.e. static expecta- 
tion values, Green's functions and the SFT grand poten- 
tial, have to be transformed back, optimization of the 
parameters is performed subsequently. 



o 
tn 
< 

£ 

E 

■ — 
as 

Q. 

a> 

"D 



0.1 



0.08 



0.06 



0.04 



0.02 



hole doped 


electron doped 




A ^ca ; 


\ 

\ 


/ /~\ VCA + bath sites " 
// 












I \ C-DMFT ; 

I i i i i I lUobd in 



0.7 0.8 0.9 



1.1 1.2 1.3 
electron density n 



FIG. 2: Superconducting (d-wave) order parameter Asc (see 
Eq. (|10[) ) as a function of the electron density n within conven- 
tional VCA (blue line, variational parameters e and hsc) and 
VCA including bath sites (red line, variational parameters e, 
Eb, V and /ib.sc)- In the case of VCA with bath sites a second 
solution with higher ground-state energy is found (dashed red 
line). Calculations have been performed for system sizes up 
to L — 6400 sites. C-DMFT results from Kancharla et al. 
(green symbols, Ref. [2(| Fig. 2) are shown for comparison. 



As an advantage of using the particle-hole transforma- 
tion, only minor changes of the standard VCA code are 
necessary. The evaluation of the self-energy functional 
can still be done using the Q-matrix technique introduced 
in Ref. |l8|, for example. 



IV. RESULTS AND DISCUSSION 

For the discussion of the results of our VCA calcula- 
tions we focus on the superconducting order parameter. 
Besides integral quantities like the ground-state energy, 
this is the observable which is sensitive to the relevant 
low-energy scale and which shows the largest discrepancy 
when comparing (plain) VCA and C-DMFT. 

To be consistent with the C-DMFT calculations of 
Kancharla et al^S we take 



\(dl c rt)\ 



(10) 



as the definition for the superconducting order parame- 
ter. The indices i and j refer to nearest-neighbor sites. 
Since translational symmetry is broken by the cluster ap- 
proach, (cQCjf) is computed for lattice sites which, in the 
reference system, would belong to the same cluster. VCA 
calculations using the reference systems given in Fig. Q] 
have been performed for a square lattice consisting of up 
to L = 6400 sites for nearest-neighbor hopping set to 
t n n — — 1) second-nearest neighbor hopping t nnn = 0.3 
and {7 = 8. 
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FIG. 3: Ground-state energy as a function of the electron 
density n within conventional VCA (blue line) and VCA in- 
cluding bath sites (red line). Dashed red line: second solution 
with higher ground-state energy. 



The results are summarized in Fig. [5] which shows the 
order parameter Age a s a function of the electron den- 
sity n for zero temperature. Apart from the normal state, 
we have allowed for a pure superconducting phase only 
and suppressed a possible antiferromagnetic phase which 
shows up close to half-filling and has been studied in 
detail in previous work^ ~ 19 i 24 In the case of the conven- 
tional VCA (blue line, reference system Fig. [TJ left), a 
Mott insulating solution is found at half-filling. With in- 
creasing doping, Age increases and reveals its maximum 
value at n ss 0.7 for hole doping and n ss 1.2 for electron 
doping. Unfortunately, we could not trace the solution 
for dopings much larger than the optimal dopings. Quite 
generally, a possible discontinuous change of the ground 
state of the finite reference cluster may result in a dis- 
continuous change of the self-energy and thus of the SFT 
grand potential. Thereby, a solution can cease to exist. 

Using VCA with bath sites (red line, reference system 
Fig. [I] bottom) the order parameter is small but remains 
finite at half-filling. The maxima of Age are found at 
n w 0.92 and n w 1.08, i.e. at significantly lower dopings 
as compared to the results of the conventional VCA. This 
comes close to the C-DMFT results of Kancharla et ali2°- 
The stationary point could not be traced for hole doping 
larger than 1 — n 0.2 and electron doping larger than 
n — 1 fx 0.3. On the hole-doped side, a second symmetry- 
broken solution could be found which comes very close to 
the C-DMFT data and shows up a second-order critical 
point with A$c ~~ * and with the symmetry-breaking 
Weiss field ft,b,SC -> at n w 0.75. This solution is found 
to exist up to n < 0.96. 

In the entire hole-doping range where two supercon- 
ducting solutions can be found, however, the one with 
the larger order parameter is more stable, i.e. for T = 



has a lower ground-state energy at the same filling. This 
is what could be expected from a (cluster) mean-field 
theory. The ground-state energies of the respective solu- 
tions are shown in Fig.[3J Note that the energy difference 
between the two solutions obtained from VCA with bath 
sites is rather small compared to the energy difference be- 
tween the results from conventional VCA and VCA with 
bath sites. It is also worth mentioning that the conven- 
tional VCA ground-state energy is higher although the 
corresponding solution has the larger order parameter. 
Both facts show that including bath sites considerably 
improves the approximation; bath sites couple to the pla- 
quette of correletad sites and yield a significant binding- 
energy gain. This gain is finite but small at half-filling 
and increases with increasing doping. 

Conceptually, the inclusion of bath sites improves the 
variational ansatz for the trial self-energy. Bath sites 
mimic the residual lattice into which the cluster is em- 
bedded by taking into account processes involving the 
sites of the cluster environment on a mean-field level. 
Physically, this means to improve the description of local 
(intra-cluster) quantum fluctuations. These additional 
fluctuations are expected and in fact seen to decrease 
the order parameter. An in principle optimal treatment 
of local fluctuations is provided by (cluster) DMFT where 
a continuum of bath sites is considered. By comparing 
the VCA with and without bath sites with the C-DMFT 
and the respective results for the order parameter (see 
Fig. [5]), we find that the essential step is already done 
by attaching a single bath site to each of the correlated 
sites. The convergence with the number of bath sites has 
already been recognized to be extremely fast in differ- 
ent contexts, see Refs. l25l - |27l . The inclusion of more and 
more bath sites must brigde the remaining discrepancies 
with the C-DMFT results since conceptually the SFT re- 
covers C-DMFT in this limit. Note, however, that for 
the practical C-DMFT calculations performed at T = 
using the Lanczos technique, a plaquette geometry with 
2 additional bath sites per correlated site has been used 
only^i 

For the conventional VCA without bath sites, the nor- 
mal state at half-filling n = 1 is described by a half-filled 
reference cluster, i.e. ml = 1. Upon doping and provided 
that the solution can be traced continuously, i.e. that 
there is no discontinuous change of the cluster ground 
state, the cluster remains at half- filling, n' = 1. This 
yields an obviously unphysical description of strongly 
doped phases. As has already been demonstrated for the 
one-dimensional Hubbard model,- 6 bath sites can help to 
overcome this problem since they serve as charge reser- 
voirs: While the reference cluster remains half-filled, the 
average occupation number n' on the correlated sites in 
the cluster is close to the electron density n of the lat- 
tice model. For a superconducting state the situation is 
somewhat different since there is no particle-number con- 
servation. Still there is a similar problem for plain VCA 
calculations without bath sites as can be seen from Fig. 
SI In the hole-doped case at n = 0.7, for example, we 
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FIG. 4: Difference \n — n'\ of the electron density in the origi- 
nal and the reference system (correlated sites only) as a func- 
tion of the electron density n. VCA calculations with and 
without bath sites corresponding to the calculations shown in 
Fig. [21 n — n' changes sign at n = 1 for the conventional VCA 
calculation (blue line) and is positive for electron doping. 



find n' » 0.88 (this corresponds to a sizeable difference 
\n — n'\ w 0.18, see figure). Adding a single bath site 
per correlated site, yields a strongly improved although 
not perfect description with a smaller difference \n — n'\ 
as shown in the figure. Conceptually, n = n' can be 
achieved for a continuum of bath sites only, i.e. for full 
C-DMFT. 

As can be seen in Fig. 2] there is a small but finite dif- 
ference | n — n'\ even at half- filling opposed to the results 
of the conventional VCA. This might correspond to the 
finite but small value of the order parameter at half- filling 
(see Fig. [5]) which then would have to be considered as an 
artifact. On the other hand, a superconducting state at 
half-filling is not unexpected for frustrated lattice models 
(i.e. for finite t nnn ) if, at the critical point for the super- 
conducting instability, the system is metallic. In fact, a 
d-wave superconducting state at n = 1 has been found, 
for example, even within plain VCA in Ref. [l9| and within 
variational Monte Carlo^ 

The translational symmetries of the original lattice are 
artificially broken by any (real-space) cluster approxima- 
tion. Since bath sites are expected to provide an im- 
proved embedding of the cluster in the residual lattice, 
however, one may expect that they help to restore the 
symmetries of the original state in a quantum-cluster ap- 
proach. It is clear that even with a full C-DMFT cal- 
culation this can be achieved only partially unless one 
uses additional (ad hoc) symmetrization schemes sub- 
sequently. Nevertheless, there is a considerable improve- 
ment as shown by our results. In Eq. (|10[) we have defined 
the order parameter Asc on neighboring sites in the orig- 
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FIG. 5: Relative difference <5a according to Eq. (|11[) as func- 
tion of n for VCA calculations with and without bath sites. 
We find the bath sites to partially compensate for the artificial 
breaking of translational symmetry. 



inal lattice belonging to the same reference cluster. Asc 
is independent of the nearest-neighbor pair for a 2 x 2 
cluster. To compare with the order parameter Asc, inter 
for neighboring sites belonging to different clusters, we 
define the asymmetry parameter 



6a 



^sc 



iSC, inter 



A 



sc 



A SC ,i: 



(11) 



Fig. [5] shows Sa as a function of the electron density 
7i. It is remarkable that the asymmetry can be reduced, 
depending on the doping, by more than one third. We 
expect that for an even stronger reduction a larger clus- 
ter would be much more efficient than adding more bath 
sites. 



V. CONCLUSION 

Several cluster-embedding approaches can be formu- 
lated within the framework of the self-energy-functional 
theory which differ with respect to bath degrees of free- 
dom. Cellular DMFT maps the lattice problem onto a 
small cluster of correlated sites to each of which a con- 
tinuous set of bath sites is attached. This ensures an 
optimal description of the local fluctuations. To access 
the zero-temperature phase diagram for Hubbard-type 
models of strongly correlated electron systems, there is 
basically the Lanczos technique only which may serve as 
a "solver" for the effective impurity (cluster) problem. 
This implies that in practice a limited number of bath 
sites can be taken into account. The plain variational 
cluster approach (VCA) employs a reference system with- 
out any bath sites at all. This approximation represents 
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a valuable counterpart to the C-DMFT which, however, 
provides an exact treatment of the local fluctuations in 
the limit of an infinitely large reference cluster only. 

An extended VCA which employs a reference system 
with a single bath site per correlated site suggests itself 
as a compromise and has been used here to study the 
doped two-dimensional Hubbard model at zero temper- 
ature. As compared to conventional (plain) VCA, we 
found a substantial gain in binding energy when attach- 
ing the bath sites. For a pure d-wave superconducting 
state, the extended VCA yields a considerably smaller 
order parameter which comes close to the predictions of 
plaquette C-DMFT. The same holds for the optimal val- 
ues for hole and electron doping, defined as maxima of the 
filling-dependent order parameter. From the comparison 
with the C-DMFT data we conclude that, as concerns 
the improved description of local correlations, attaching 
a single bath site per correlated site does the main part 
of the job. Clearly, larger clusters exceeding a plaquette 
are desirable to improve the description of spatial corre- 
lations. 

Methodically, bath sites serve as charge reservoirs and 
in this way yield an average occupation on the correlated 
sites in the cluster which, opposed to plain VCA, is much 
closer to the band filling in the original lattice model 



although the total cluster filling still stays at half-filling in 
the normal state (this is different for the superconducting 
phase where the particle number has no definite value). 
Bath sites also help to partially restore the translational 
symmetry that is artificially broken by the (real-space) 
quantum-cluster concept. 

Physically, the comparison of the different variational 
cluster approximations shows that local temporal cor- 
relations in addition to non-local spatial correlations 
are important to understand the size and the dop- 
ing dependence of the d-wave superconducting order 
parameter. The corrections due to local correlations 
are most pronounced for the optimal-doped and the 
overdoped regime. Here, local quantum fluctuations are 
found to substantially decrease the order parameter. 
This could be attributed to a Kondo-type screening 
of local magnetic moments which in turn reduces spin 
fluctuations. 
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